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MOTIVATION 

Accurate and Efficient Prediction of Emissions 


1. Accurate Prediction of Emissions From Combustion 
Devices Requires Treatment of Finite-Rate Kinetics 

2. The Effect of Turbulent Fluctuations in Velocity, Energy, 
Composition, etc. on Finite-Rate Chemical Kinetics Must 
be Modeled 
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TURBULENCE/CHEMISTRY INTERACTIONS 
Possible Approaches 


• Neglect Fluctuations 
+ Simple 

- ignores Effect of Turbulence 

• Eddy Break Up 

+ Simple 

Assumes Fast Chemistry 

Mean Density, Temperature Must Still Be Modeled 

• Prescribed PDF 
+ Efficient 

Limited to Fast Chemistry or Single Step Reaction 

• Composition PDF 

+ Finite-Rate Multi-Step Kinetics 
Expensive 
Gradient Diffusion 

• Velocity-Composition PDF 
+ More Accurate 

More Expensive 


PARTICLE REPRESENTATION 

A Solution Method for a Large Number of 
Independent Variables 


• Computational Requirements Increases Exponentially 
With Dimensions for Finite Difference Methods 


• Computational Requirements Increase Linearly 
With Dimensions for Monte Carlo Methods 
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COMPOSITION PDF SOLUTION 

Stochastic Lagrangian Particle Simulation 


Particle Composition and Position Changed to Model Transport of 
Joint PDF 


• Mean Convection 

- Move Particles Between Cells 


* Chemical Reactions 

Lookup Table Holds Composition Change 


• Turbulent Diffusion 

- Exchange Particles Between Cells 



• Molecular Mixing 

- Particle Interaction Changes Composition 


COUPLING 

PDF Solution is Separate Module 


CFD-ACE 


Monte Carlo 
PDF 


u, v, w, k,s 


u, v, w.'p 

f (% , ...y„) 

k,e 

"p 
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CHEMICAL KINETICS 

Reduced Models are Used 
Hydrogen: 2H 2 + 0 2 <=>2H 2 O 

CO: C0 + H 2 0oC0 2 + H 2 

2H 2 + 0 2 <=>2H 2 0 

Methane: CH« + 2H + 2H 2 0 -» CO + 4H 2 

C0 + H 2 0oC0 2 + H 2 
2H 2 + 0 2 -»2H 2 0 
3H 2 +0 2 «2H 2 0 +2H 

Hydrocarbon :C n +(y-)0 2 -> n CO +(n +1) H 2 
C n »2n*2* n H 2 0 -* n CO +(2n +1) H, 

C0 + H 2 0«C02+ h 2 

2H 2 +0 2 »2H2 0 

Thermal NO: N 2 +0«N0 + N 
N +0 2 »N0 + 0 
N +OH <=> NO + H 


RESULTS TO BE PRESENTED 


• Jet Diffusion Flame (Hydrogen with Helium Dilution) 

• Bluff Body Stabilized Flame (H 2 /CO) 

• Piloted Jet Diffusion Flame (Methane) 

• Generic Gas Turbine Combustor (Propane) 


122 





HYDROGEN JET DIFFUSION FLAME 


Illustration of Experiment at 
Sandia National Lab 


COUUST 

RJ7W 





Re = 10 4 


Fuel 

100% H 2 
80% H 2 , 20% He 
60% H 2 , 40% He 


60% HYDROGEN FLAME 

Scatter Plots of Mixture Fraction and 
NO Mole Fraction 
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UNDILUTED HYDROGEN FLAME 
Conditional Averged NO Mole Fraction 






HYDROGEN DIFFUSION FLAME 
Dilution Effects on Emmisions Index 



0 5 10 15 20 25 30 35 40 

Dilution (% Ho) 
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BLUFF BODY STABILIZED DIFFUSION FLAME 
Illustration of Experiment of Correa and Gulati 


Inlet Air Flow 


27.5%CO 
32.3%H2 
40.2% N2 


15 cm x 15 cm windtunn«4 
top w«a 



bottom wad 


BLUFF BODY STABILIZED DIFFUSION FLAME 
Composition PDF Predicts Mean Values as 
well as Velocity-Composition PDF 
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PILOTED JET DIFFUSION FLAME 
illustration of Experiment of Masri et.al. 


Experimental Conditions 


Flame 

Fuel jet Velocity 

L 

41m /s 

6 ' 

48 m/s 

M | 

55 m/s 


ntnaiD 



PILOTED JET DIFFUSION FLAME 
Good Agreement with Experimental Data 





x/D = 20 


x/D = 30 
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PI LOTED JET DIFFUSION FLAME 

More Accurate Prediction with Monte Carlo PDF 

No PDF Prescribed PDF Monte Carlo PDF 



Flame B at x = 20D 


GENERIC GAS TURBINE COMBUSTOR 
Pratt & Whitney Four-Nozzle Sector 
Combustor Tested at Wright Laboratory 



67,840 Cells 
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MO NTE CARLO PDF COMBUSTOR CALCULATION 
Stochastic Particle Traces """" " 



VERTICAL PLANE THROUGH CENTER OF 

FUEL INJECTOR 

Mean CO Mass Fraction Countours 
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RUN TIME AND MEMORY 


3D Combustor Calculation 
(68,000 cells) 


Conventional CFD 


CPU Time 

20 hours 

Memory 

80 MBytes 


Monte Carlo PDF 


CPU Time 

100 hours 

Memory 

120 MBytes 


Parallel PDF (Projected) 


CPU Time 

25 hours 

25 hours 

25 hours 

25 hours 

Memory 

30 MBytes 

30 MBytes 

30 MBytes 

30 MBytes 


CPU Time for IBM RS/6000 Model 560 


CONCLUSIONS 


• Monte Carlo PDF Solution Successfully Coupled with 
Existing Finite Volume Code 

- Minor Changes to Finite-Volume Code 

- Can be Coupled with Other Codes 

• PDF Solution Method Applied to Turbulent Reacting 
Flows 

- Good Agreement with Data for 2D Case 

- Demonstration of 3D Elliptic Flow 

• PDF Methods Must be Run on Parallel Machines for 
Practical Use 
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